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Antiferromagnetic Heisenberg spin chains in a sufficiently strong magnetic field are Luttinger liq- 
uids, whose parameters depend on the actual magnetization of the chain. Here we present precise 
numerical estimates of the Luttinger liquid dressed charge Z, which determines the critical expo- 
nents, by calculating the magnetization and quadrupole operator profiles for 5 = 1/2 and S = 1 
chains using the density matrix renormalization group method. Critical amplitudes and the scat- 
tering length at the chain ends are also determined. Although both systems are Luttinger liquids 
the characteristic parameters differ considerably. 

PACS numbers: PACS numbers: 75.10.Jm, 75.40.Cx, 75.10.-b 
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I. INTRODUCTION 

The one-dimensional (ID) antiferromagnetic Heisen- 
berg chain 



N-l 



N 



(1) 



is one of the most thoroughly investigated paradigms of 
interacting many-body systems. It is well-known that 
for zero magnetic field the low-energy physics depends 
very much on the spin length 5*. As Haldane predicted 
in 1983 ?i integer 5 chains possess an energy gap A g (S) 
above the ground state, the size of which vanishes expo- 
nentially as 5 — > 00, whereas half- integer 5 chains are 
gapless. The presence of the Haldane gap in integer spin 
chains implies the occurrence of a critical field h c = A g 
beyond which the system gets magnetized. As the gap 
collapses at the critical field the ground state structure 
changes adequately, and in an already magnetized state 
there remains no conceptual difference in the low-energy 
properties between integer and half-integer chains. 

Partially magnetized antiferromagnetic Heisenberg 
chains at low energies are expected to be one-component 
Luttinger liquids 2 (LLs) irrespective of the spin length 5. 
This is known rigorously in the 5 — 1/2 case where the 
model can be analyzed using the Bethe Ansatz^ The va- 
lidity of the Luttinger liquid description for 5 > 1 chains 
and coupled spin- 1/2 ladders has been investigated and 
confirmed by analytical^*& and numerical^&S methods. 
Although there is a theoretical possibility of finding gap- 
ful behavior (magnetization plateau) at special values of 
the magnetization, such a scenario does not seem to be 
realized in pure nearest-neighbor Heisenberg chains (see 
Ref. for 5 = 3/2). Neither would we expect multi- 
component LLs, which can otherwise occur for more gen- 
eral (higher-order) couplings iii 

The LL concept bears a strong relationship to c = 1 
conformal field theory (CFT),- 12 In fact the LL is hardly 
more than the CFT adapted to the situation with two 
Fermi points 2kp momentum apart. The one-component 



LL is a three parameter theory^ In the present context 
the first parameter, the location of the Fermi points ±/cf, 
is determined by the magnetization through 



2k F = 2tt(5 - to) 



(2) 



where to = S^ ot /L is the (bulk) magnetization. The 
second parameter, the Fermi velocity v is just an en- 
ergy scale, whereas the third parameter, to be called 
the "dressed charge" Z, determines the universality class 
and the critical exponents. (In the bosonization litera- 
ture the notation K = Z 2 is standard.) The traditional 
LL parameters, the velocities for charge and current ex- 
citations can be expressed with Z as w c hargo = v/Z 2 , 

^current — ^ Z ^™ 

All three parameters are functions of the actual mag- 
netization m and the spin length 5. The LL parame- 
ters of the 5 — 1 chain in a magnetic field have already 
been determined to some extent using numerical meth- 
ods. Sakai and Takahashi diagonalized small finite chains 
up to N = 16 with periodic boundary condition and used 
the prediction of CFT on finite size energy spectra to esti- 
mate the critical exponents £ A similar method was used 
by Usami and SugaS for the 5=1/2 ladder which be- 
haves as a Haldane-gap system for strong ferromagnetic 
interchain interaction. More recently Campos Venuti et 
al. used the DMRG to compute directly the transverse 
two-point function G xx (r) — (S x (0)S X (r)) on chains with 
N = 80 and fitted these data to the expected asymptotic 
form of G XX H 

In this paper we apply an alternative method. We cal- 
culate precise numerical estimates to the critical expo- 
nents and by this to the LL dressed charge Z via a direct 
numerical determination of the magnetization profiles in 
finite open chains. The magnetization profile is, by defi- 
nition, the positional dependence of the local magnetiza- 
tion m n — (5^) in an open chain with some well-defined 
boundary condition. The form of the profile depends on 
the applied boundary condition, and may involve sur- 
face exponents as well whenever the boundary condition 
on the left and right ends differ. In order to simplify 
the expected behavior we will only consider cases where 



2 



the boundary condition is identically open (free) on both 
ends; in this case only bulk exponents come into play^ In 
a semi-infinite chain, far from the chain's end, the mag- 
netization profile is expected to decay to its bulk value 
m algebraically as 



m(r) ~ m + A m cos(2fc F r + (j)) r r '^ 2 1 



(3) 



where A m is a non-universal amplitude, is a phase 
shift, rj z is the (bulk) critical exponent, defined through 
the translation invariant longitudinal two-point function 
G zz (r) = (S z (0)S z {r))~m 2 ~ A cos(2kpr + (f>) r -17 * , and 
kp is determined as a function of m by Eq. J2}. The form 
of the Friedel oscillation in Eq. (|3} contains indirect in- 
formation about the LL parameter Z ', and this can be 
exploited in a numerical procedure to calculate precise 
estimates. This program has been carried out for the 
5 = 1/2 chain and ladders in Ref. |3, for the Kondo lat- 
tice model in Ref. Il6l for the Hubbard model in Ref. 0, 
and for the t — J model in Ref. 

While the algebraic decay with rj z /2 is a standard con- 
sequence of criticalityji^ the 2kF oscillation is a special 
LL feature, which stems from the fact that we work with 
two families of (chiral) CFT operators living around the 
two Fermi points in fc-space. Without the 2kp oscillation 
it is a standard exercise in CFT to derive the shape of 
the magnetization profile in a strip geometry of width L 
by applying the logarithmic mapping. 15,2 i This transfor- 
mation, together with the proper account for the 2kp LL 
term yields a prediction for the magnetization profile of 
a finite Luttinger liquid segment of length L with open 
boundary condition on both ends&iS 
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(4) 

Note that there is no need to introduce explicitly a phase 
shift 4> as in the semi-infinite case Eq. ©, since the sym- 
metry of the profile with respect to r — > L — r implies 
that = —kpL. 

The predicted LL profile in Eq. Q is based on the 
continuum limit. However, for finite lattice systems such 
as the Heisenberg chain, we expect corrections. Confor- 
mal invariance and thus Eq. (0} is only expected to be 
valid asymptotically in the large L limit with r satisfying 
<< r << L. Clearly, in a strict sense the magneti- 
zation cannot diverge at the chain ends since m n < 5. 
Phenomenologically this natural cut-off at the boundary 
acts as an effective impurity put into the CFT model 
at the system ends. This defines an effective "scatter- 
ing length", SN = 0(1), associated with the bound- 
ary, suggesting that we should replace the naive system 
size L — N in Eq. 0} with an "effective" system size 
L — > N + 2 SN. Based on this intuitive argument in the 
following we will work with a slightly modified Ansatz 
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N + 2SN , 
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(5) 



where 5N is a free fitting parameter. When SN = Eq. 
(JSJ) is simply the discrete version of Eq. 10} (note that 
now the system is defined for n — 1, . . . , N). 

In the lack of a consistent scheme to calculate correc- 
tions, the above amendment to the fitting formula is not 
rigorous. Justification can stem from a direct compari- 
son with exact results. In the following we first discuss 
in detail the exactly solvable 5=1/2 chain as a bench- 
mark case. We demonstrate that the Ansatz in Eq. J5J is 
capable of reproducing the (numerically) exact 5=1/2 
critical exponents to high precision even with SN = 0. 
When SN is also fitter for, the accuracy improves by an 
order of magnitude. 



II. THE 5=1/2 CHAIN: BETHE ANSATZ 
RESULTS 

The 5=1/2 Heisenberg chain in a magnetic field can 
be solved exactly by the Bethe Ansatz (BA) method. 2 - 
The BA itself cannot yield the correlation functions and 
critical exponents directly, but assuming conformal in- 
variance the operator content can be read off from the 
structure of the low-energy excitations above the ground 
stated These latter can be determined by a systematic 
calculation of the finite size corrections to the L = oo 
BA equations (assuming periodic boundary condition). 22 
The energy and momentum of the lowest energy (pri- 
mary) states (with respect to those of the ground state, 
Eg and P g ) can be cast into the form 



2ttv 



SE = E & -E g = — (A+ + A-), 
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5P = p a - Pg = Q + — (A+-A-), 



(6) 
(7) 



where a = {d, 1} is a shorthand for two integer topo- 
logical quantum numbers labeling the states. The 0(1) 
momentum term is 



Q = Qa = 2kpd + irl, 



(8) 



with kp defined in Eq. @. The "conformal dimensions" 
A ± read 
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(9) 



where Z is the dressed charge. The topological quan- 
tum numbers have a direct physical interpretation in the 
LL representation: I (d) denotes the number of fcrmions 
added to (transferred from the left Fermi point to the 
right in) the band. 

The critical exponents appearing in the one- and two- 
point correlation functions can be expressed with the con- 
formal dimensions A ± . In general a physical operator 
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FIG. 1: The dressed charge as a function of the magnetization 
in an S — 1/2 chain; o denotes results using the naive fit 
SN = 0, + denotes fits where <5iV is optimized. Solid line is 
the Bethe Ansatz result. Inset shows the critical exponents 
and the numerically determined amplitude A m . At m = 
the dressed charge is Z = l/y2. 



such as S z decomposes into all operators which are not 
forbidden by conservation (selection) rules. In particu- 
lar, for S z the number of fermions in the band is con- 
served, thus necessarily I — 0. The asymptotic decay is 
determined by the operator which has the lowest critical 
exponent, i.e., d — 1 giving 



r lz = 2(A+ + A-)\ l ^ d=1 = 2Z*, 



(10) 



and an oscillation 2k p by Eq. (JHj). The transverse critical 
exponents r) x is determined by the operator I = 1, d = 0, 
i.e., 



2(A H 
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(11) 



with the associated oscillation frequency Q = ir. Equa- 
tions lfffi|) and ifTT|) imply rj z r] x = 1. 

In the BA approach the dressed charge Z = Z(B) is 
determined by a set of integral equations for the den- 
sity of rapidities p{x), and the dressed charge function 



where 



p(x) 
Z(x) 



K(x — x )p(x )dx , 



K(x-x')Z(x')dx', 
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(12) 
(13) 

(14) 



The limits of integration are defined implicitly through 
the constraint 



p(x')dx = 1/2 — m. 



(15) 



FIG. 2: Magnetization profiles determined by the DMRG for 
the S = 1/2 chain with N = 120, M = 160 and magnetization 
m = S z ot /N w 1/3. Only half of the chain shown - the other 
half is mirror symmetric. 



The coupled integral equations in Eqs. I|12|) . (|13|) and 
(|15fl are to be solved numerically^ First, assuming B is 
given, Eq. (|12|) is solved for p. This is inserted into Eq. 
(|15fl to find the associated value of m. When m is given, 
as in our case, this procedure can be iterated to find B as 
a function of m at arbitrary precision. Finally Eq. (|15fl is 
solved for the function Z(x) and the dressed charge Z = 
Z(B) is determined. The numerically determined dressed 
charge as a function of the magnetization is plotted in 
Fig. ^ I n the next section this will serve as a reference 
curve to check the accuracy of the fitting Ansatz in Eq. 



III. THE 5=1/2 CHAIN: DMRG RESULTS 

In order to test the fitting formula Eq. (JSJ) we calcu- 
lated the magnetization profile in the ground state of 
finite chain segments with open boundary condition us- 
ing the density matrix renormalization group method»2s. 
The chain length was N — 120 and we kept M = 160 
states. The truncation error was found to be in the range 
1 - P M = 10~ 13 - KT U . We applied the finite lattice 
algorithm with four iteration cycles which was found suf- 
ficient for convergence. An example of the magnetization 
profile as determined by the DMRG is depicted in Fig. 
® 

It is worth discussing the fitting procedure itself in 
some detail. The Ansatz in Eq. JSJ has five parame- 
ters: m, A m , fci?, 5N and rj z . Although S* oi is a con- 
served quantity, and thus can be set to a given value 
in the DMRG, due to the open boundary condition the 
control over the exact (bulk) value of m is lost. While 
it is true that for long enough chains m will be close to 
S* ot /N, the finite size deviation should be tracked during 
the fitting procedure. Similarly, although kp is a well- 
specified function of m in the bulk [see Eq. @], it was 
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FIG. 3: Stability of local fit parameters as a function of local 
position in the chain for S — 1/2, L — 120, m ~ 1/4. Local 
fit parameters at position n axe defined by fitting for sites 
i — n, . . . , n+3A, with A = 2iv/kF = l/(S—m) the wavelength 
of the oscillations. The highest stability of the parameters is 
achieved for SN « 0.4. Equation is satisfied up to 2.10 -4 . 

found advantageous to keep it as an independent fit pa- 
rameter, and only use Eq. a posteriori as a consistency 
check. On the other hand, it is better to avoid fitting on 
SN directly. Instead the best working alternative seems 
to be making a four-parameter fit on to, A m , kp and ?y 2 , 
while SN kept fixed. The optimal value of SN is the one 
which yields the highest stability with respect to local 
fits, i.e., calculating the four fitting parameters from a 
small number of sites at different locations in the chain. 
An example of this procedure is shown in Fig. [21 We 
found that the optimal value of SN is a weak function 
of m, being about SN s» 0.5 for to = 0, and decreasing 
monotonically to SN as 0.4 for to = 1/2. This is more 
or less consistent with the value SN = 1/2 used in the 
bosonization approach of Ref. 0- 

Having obtained r\ z using the above fitting procedure, 
the fundamental quantity of the theory, the dressed 
charge Z, can be calculated from Eq. I|1U|) . Figure Q] 
shows the numerically determined dressed charge as a 
function of to. The relative error of the fitting proce- 
dure is under 0.1%, except for very small to values where 
logarithmic corrections to the fitting formula and Z are 
expected, and the fitting procedure loses stability (see 
the error bar in the figure at to = 0.03). Otherwise, the 
agreement with the exact values is very good and stable 
for 0.03 < to < 1/2. For comparison Fig.njalso shows the 
estimate of Z when the naive scattering length SN = 
is used. Even with this the error is within 1%, except for 
very small to. 

We conclude that for S = 1/2 the LL Ansatz for the 
magnetization profile, Eq. JSJ, is a very efficient tool in 
calculating the critical exponents (dressed charge). The 
final result is highly accurate already with SN = 0, but 
an additional increase in precision can be achieved by 
adding the scattering length as an extra fit parameter. 



IV. THE S = 1 CHAIN 

Encouraged by the success of the fitting procedure for 
5 = 1/2, in this section we apply it to the S = 1 Heisen- 
berg chain. The forthcoming analysis is not intended 
to prove that the S = 1 chain in its magnetized regime 
is a Luttinger liquid - this has been done convincingly 
alreadyi 4 i 5 i 6 i 7 i 14 Instead, our starting point is the assump- 
tion that we have a one-component Luttinger liquid, and 
then use LL theory and the numerically calculated mag- 
netization profiles to determine the value of Z and the 
critical exponents with high accuracy. 

Although the spin-1 problem cannot be treated rigor- 
ously, there are two limits where clear theoretical predic- 
tions exist. In the high magnetization limit near satu- 
ration the physics can be understood by regarding the 
system as a dilute gas of magnons created in the fer- 
romagnetic vacuum^ Magnons behave as bosons with 
short-range repulsive interaction. In the dilute limit the 
exact form of the interaction is irrelevant and an essen- 
tially hard-core boson description becomes valid. In ID 
hard-core bosons are equivalent to free spinless fermions, 
which imply an LL parameter Z = 1 with correlation 
exponents r\ z = 2 and rj x = 1/2 as m — ► 1. 

There is a similar argument in the to — ► limit. Be- 
low the critical magnetic field h c < A g , where A g is 
the Haldane gap, the elementary excitations are mas- 
sive spin-1 bosons^ At h c these bosons condensate, but 
due to the interboson repulsion the magnetization only 
increases gradually. Near to > the boson density is 
low and by the above argument we again obtain Z = 1, 
rj z = 2 and ij x = 1/2. Note that this argument only works 
for S = 1 with a Haldane gap. For S =1/2 for which the 
nonmagnetized ground state is massless, the elementary 
excitations form a strongly interacting dense gas with 
Z = 1 / V2 as is given by the Bethe Ansatz, and the SU(2) 
symmetry at m = 0. Alternative theories for S = 1, 
based on either a Majorana fermion representation* or 
on the bosonization of spin- 1/2 ladders^ lead to a similar 
conclusion in this limit. 

Although the magnetized S — 1 chain is a Luttinger 
liquid, this classification only applies at low energies and 
long distances. Indeed, at higher energies and shorter dis- 
tances the S = 1 chain produces features which cannot 
be understood within the framework of LL theory. These 
features, absent in S = 1/2 chains, stem from the addi- 
tional degrees of freedom staying massive for S > 1/2. 
These degrees of freedom have signature both in the en- 
ergy spectrum and correlations, and manifest themselves 
in the numerical, finite chain calculations. Their ori- 
gin can be easily understood in the low-magnetization 
limit. At zero magnetic field the system possesses an 
energy gap, the Haldane gap. The lowest excited states 
form a triplet branch with a minimum energy at momen- 
tum k = tt. The operator Sf has large matrix elements 
between the ground state and the S* ot = component 
of this triplet. This leads to an exponentially decaying 
alternating (antiferromagnetic) behavior in the longitu- 
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FIG. 4: Magnetization profiles determined by the DMRG for 
the 5" = 1 chain with TV = 120, M = 160 at the bulk mag- 
netization m ~ 0.097. Only half of the chain shown - the 
other half is mirror symmetric. The LL oscillation is clearly 
visible above n > 20. Solid line is the LL fit. Inset shows the 
end of the chain on a log-lin scale featuring the exponentially 
decaying initial oscillation which crosses over to LL behavior 
at n w 12. 
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FIG. 5: Dressed charge Z determined numerically for the 
S — 1 chain. Open symbols denote results for L — 120, closed 
symbols for L — 240. The estimated error is comparable to 
the symbol size. Solid line is only a guide to the eye. Inset 
shows the longitudinal and transverse critical exponents r} z 
and r) x , resp. 



dinal correlation function. When the magnetic field is 
switched on, the Zeeman energy splits the triplet branch, 
and at h = A g the 5 t z ot = 1 component at k — tt crosses 
over with the ground state. However the 5 t z ot = com- 
ponent remains in the spectrum (at energy ss A s ) and 
still contributes to the short-range longitudinal correla- 
tion functions. As a consequence the two-point function 
shows a crossover from a seemingly exponential decay on 
short distances caused by the massive mode to an alge- 
braic decay determined by the soft, LL mode on longer 
distances. 

There is a similar effect in the one-point function we 
consider here. Near the chain's end there is an exponen- 
tially localized effective 5=1/2 degree of freedom, the 
so-called " end spin" , which also survives in the magne- 
tized regime, at least when the magnetization is not too 
high. 26 Its presence produces a crossover from exponen- 
tial decay to algebraic decay in the one-point correlation 
function as is illustrated in Fig. As the magnetization 
increases the massive modes rise in energy, and have a less 
and less significant impact on the low-energy physics. At 
the same time the end spins gradually dissolve and dis- 
appear in the bulk as was observed by Yamamoto and 
MiyashitaiS 6 . 

In order to measure the critical exponent rj z we used 
the DMRG algorithm with N = 120 and M = 160 with 
five iteration cycles. The truncation error varied in the 
range 1 — P m = 10~ 12 — 10~ 8 , the calculation being more 
precise in the high magnetization regime. A limited num- 
ber of runs with N = 240, M — 300 was done to check the 
numerical precision and to obtain results at points where 
longer systems were needed. We applied the fitting pro- 
cedure described above and care was taken to stay in the 
bulk of the chain sufficiently far from the ends. Since 



for to = the correlation (localization) length of the end 
spin is about 6 lattice sites, which becomes even shorter 
for m > 0, a chain length N = 120 was found sufficient. 

We observed that the fitting procedure is somewhat 
less accurate here than for 5 = 1/2, meaning that finite 
size corrections to the Luttinger liquid profile Eq. (JjjJ 
are more important in for S = 1. The fitting procedure 
becomes especially unreliable below to « 0.05 and around 
to « 0.5. In the former case the increasing wavelength, 
which becomes comparable to the system size, whereas in 
the latter case the vanishing prefactor cos[2fci?(n — (N + 
l)/2)] — > can be blamed for the numerical difficulty. 
Far enough from these problematic regions the relative 
error of the calculated exponent was estimated to be less 
than 1%. 

Beyond measuring the local magnetization profile m„, 
for S — I there is another, independent quantity whose 
profile can also be measured easily. This is the local 
quadrupole moment q n = ((S*) 2 ). Note that for 5=1/2 
this quantity is redundant and thus carries no additional 
information. For 5=1, however, the quadrupole profile 
provides us an alternative way to measure the critical 
exponent which in many cases became even more precise 
than the one obtained through the magnetization profile. 
The quadrupole profile is expected to behave according 
to the same scaling form Eq. (J5J with the replacement 
to -> q, A m -> Ag. 

Figure [S] shows the dressed charge determined using 
Eq. I|10fl from the measured rj z exponent of magnetiza- 
tion and quadrupole profiles. The inset also shows the 
corresponding critical exponents r\ z and r\ x as a function 
of the bulk magnetization. We see that for m — > 1 the 
predicted value Z = 1 is reached very rapidly. We an- 
alyzed the m-dependence close to m = 1 and found it 
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FIG. 6: Magnetization and quadrupole moment amplitudes 
for the 5=1 chain. Open symbols denote results for L = 120, 
closed symbols for L = 240. Inset shows the best estimate of 
8N. 

to be describable with a power law with a rather large 
exponent around 4, although a scaling faster than any 
power law cannot be excluded either. Due to the lack of 
sufficient numerical precision near m — 1 we were unable 
to resolve this question reliably. 

For m — ► the conclusion is also somewhat vague be- 
cause the m-dependence is very steep and due to nu- 
merical difficulties we were unable to approach this limit 
closely enough. However, a Z = 1 value predicted by 
the theory seems highly consistent with our data. As- 
suming this, our numerical values seem to indicate an 
m-dependence Z — 1 ~ m a , with a = 0.5 ± 0.1. Knowing 
that to scales above the gap asm~ yjh — A g , 27 this 
would imply Z — 1 ~ (h — A g )P for small h with j3 w 1/4. 

Between the two limits Z is larger, reaching its largest 
value Z w 1.21 at m ~ 0.36. This behavior should be 
compared to that of 5 = 1/2, where Z < 1 in the whole 
regime (see Fig. H). Earlier data on the critical expo- 
nents available in the literature, which were determined 
by other methods such as the finite-size scaling analy- 
sis of the energy spectrum^ or the staggered structure 
factor^ are consistent with our results. 

The numerically determined critical amplitudes are de- 
picted in Fig. The amplitude of the magnetic moment 
fluctuations A m increases monotonically as a function of 
to. However, the rate of increase is not smooth as it is 
seen in the figure. A m changes sign somewhere close to 
to = 1/2. Where A m w the observable fluctuations are 
governed by the next smallest critical exponent, and thus 
a precise measurement is beyond our method. In contrast 



with this, the amplitude of the quadrupole fluctuations 
A q remains positive in the whole magnetized regime. We 
observe that it decreases for small to, reaching its min- 
imal (still positive) value at to ps 0.36. Above this it 
increases and saturates for m — ► 1. 

Finally, it is interesting to note that the optimal value 
of the scattering length SN changes considerably as m 
varies. As the inset of Fig.[f)]shows that 5N is around 0.5 
for small to, then decreasing to 5N w —0.5 at to w 0.7, 
then increasing again to 5N « at m = 1. There is 
a relatively large uncertainty in the optimal value de- 
termined, but this imprecision has little impact on the 
estimated value of Z and the correlation exponents. 

V. SUMMARY 

In this paper we have analyzed the critical fluctuations 
in open 5 = 1/2 and S = 1 Heisenberg chains in their 
magnetized regime. The low energy physics is a one- 
component Luttinger liquid in both cases. We have calcu- 
lated the LL's characteristic dressed charge and the am- 
plitude of the leading critical fluctuations. Our method 
consisted of determining numerically the magnetization 
and quadrupole operator profiles and applying a fitting 
procedure based on conformal invariance. The method 
has been thoroughly tested on the Bethe Ansatz solvable 
5 = 1/2 chain, confirming its reliability and high pre- 
cision. In the 5=1 case, where the exact solution is 
unknown the method provided high-precision estimates 
of the critical exponents, justifying and complementing 
earlier results derived by alternative methods. We also 
determined critical amplitudes which has much less been 
studied so far for these systems. 

Beyond calculating the characteristic parameters with 
high accuracy our results also allow us to make a detailed 
comparison between the 5 = 1/2 and 5=1 chains. Al- 
though both are Luttinger liquids, the respective dressed 
charges as functions of the magnetization differ quite con- 
siderably. The only limit where the two systems become 
equivalent is the full saturation limit, m — > 5, where 
Z — ► 1. Otherwise, for 5 = 1/2 the dressed charge is 
Z < 1, while for 5 = 1 it is Z > 1. There are interesting 
differences in the behavior of the critical amplitudes and 
the scattering lengths at the chain ends, as well. 

We thank F. Igloi and J. Solyom for valuable discus- 
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